\chapter{Kurzschlussversuch}
Beim Kurzschlussversuch wird die zu vermessende Maschine von einer zweiten Maschine angetrieben. Aus diesem Grund besteht zwischen beiden Maschinen eine mechanische Kupplung. Die Anschlussklemmen der zu vermessenden Maschine sind kurzgeschlossen und es wird jeweils der Strom in jedem Strang, die Drehzahl und das Drehmoment zwischen angetriebener und antreibender Maschine gemessen.
\section{Theorie}
\label{kurz_theo}
\begin{figure}[htb]
	\centering
  \includegraphics[width=0.8\textwidth]{auswertung_kurzschluss/ESB_kurz.png}
	\caption{Ersatzschaltbild je Phase einer PM-Synchronmaschine für den Kurzschlussversuch}
	\label{Ersatzschaltbild_kurz}
\end{figure}
Für den Kurzschlussversuch ergibt sich das Ersatzschaltbild je Phase nach Abbildung \ref{Ersatzschaltbild_kurz}. Die Polradspannung $\underline{U}_\+{p}$ sorgt dafür, dass ein -Kurzschlussstrom $\underline{I}_\+{s,k}$ fließt, welcher sich direkt aus dem Spannungsumlauf ergibt:
\begin{equation}
 \underline{I}_\+{s,k}= -\frac{\underline{U}_\+{p}}{R_\+{s}+j \omega_\+{s} (L_\+{h}+L_\+{s,$\sigma$})}
\label{eq:ku1}
\end{equation}
Für den Betrag des Kurzschlussstroms $I_\+{s,k}$ gilt damit:
\begin{equation}
 I_\+{s,k}=|\underline{I}_\+{s,k}|= \frac{U_\+{p}}{\sqrt{R_\+{s}^\+{2}+\omega_\+{s}^\+{2}(L_\+{h}+L_\+{s,$\sigma$})^\+{2}}}
\label{eq:ku2}
\end{equation}
Die Statorkreisfrequenz $\omega_\+{s}$ steigt linear mit der Drehzahl an. Bei hohen Drehzahlen ist somit der Einfluss des Strangwiderstandes $R_\+{s}$ auf den Kurzschlussstrom vernachlässigbar. Der Betrag und die Phasenlage des Kurzschlussstroms werden daher hauptsächlich durch die Polradspannung, sowie die Induktivitäten $L_\+{h}$ und $L_\+{s,$\sigma$}$ bestimmt. Abbildung \ref{Zeigerdiagramm_kurz} stellt im Kurzschlussfall das Zeigerdiagramm je Phase dar. Es ist zu erkennen, dass der resultierende Stromzeiger im Wesentlichen in der d-Achse liegt. Damit lässt sich durch Umstellen von Gleichung \ref{eq:ku2} und Einsetzen der, aus dem Leerlaufversuch bekannten und zur Drehzahl proportionalen, Polradspannung die Induktivität $L_\+{d}$ bestimmen:
\begin{equation}
L_\+{d}=L_\+{dh}+L_\+{s,$\sigma$}=\frac{\sqrt{U_\+{p}^\+{2}-R_\+{s}^\+{2}I_\+{s,k}^\+{2}}}{\omega_\+{s}I_\+{s,k}}
\label{eq:ku3}
\end{equation}
\input{auswertung_kurzschluss/zeiger}
%\begin{figure}[htb]
	%\centering
  %\includegraphics[width=0.5\textwidth]{auswertung_kurzschluss/Kurzschluss_Bild.png}
	%\caption{Zeigerdiagramm Kurzschlussversuch}
	%\label{Zeigerdiagramm_kurz}
%\end{figure}
Das von der angetriebenen Maschine abgegebene elektrische Drehmoment $M_\+{e,k}$ kann über die Luftspaltleistung berechnet werden. Die Luftspaltleistung setzt sich aus der Polradspannung und dem Anteil des Kurzschlussstroms, der die gleiche Phasenlage aufweist, zusammen. Da die Polradspannung in der q-Achse liegt, ist dies der Imaginärteil. Die Zusammenhänge ergeben sich wie folgt: 
\begin{equation}
M_\+{e,k}=\frac{P_\+{$\delta$}}{2 \pi n}=\frac{3 U_\+{p}}{2 \pi n}\operatorname{Im} \{ \underline{I}_\+{s,k} \} = \frac{3 U_\+{p}^\+{2}}{2 \pi n}\frac{R_\+{s}}{R_\+{s}^\+{2}+\omega_\+{s}^\+{2}(L_\+{h}+L_\+{s,$\sigma$})^\+{2}}
\label{eq:ku4}
\end{equation}
Gleichung \ref{eq:ku4} kann genutzt werden, um die Abhängigkeit zwischen der Drehzahl und dem Kurzschlussmoment $M_\+{k}$, welches an der Drehmomentmesswelle zwischen den beiden Maschinen gemessen wird, zu bestimmen. Dabei ist zu beachten, dass die Statorkreisfrequenz und die Polradspannung linear mit der Drehzahl ansteigen. Wird beispielsweise der, im Leerlaufversuch bestimmte, Faktor $k_\+{Up}$ eingesetzt, so ergibt sich nach Kürzen mit der Drehzahl in \si{\per\minute} der folgende Zusammenhang:
\begin{equation}
M_\+{e,k}= \frac{90 k_\+{Up}^\+{2}n}{1000^\+{2}\pi}\frac{R_\+{s}}{R_\+{s}^\+{2}+2\pi(\frac{n p}{60})^\+{2}(L_\+{h}+L_\+{s,$\sigma$})^\+{2}}
\label{eq:ku5}
\end{equation}
Bei geringen Drehzahlen dominiert der Zähler, sodass das Kurzschlussmoment zunächst nahezu linear über der Drehzahl ansteigt. Bei hohen Drehzahlen nimmt das Kurzschlussmoment etwa hyperbolisch ab, da nun der Nenner dominant ist. Die gemessene Kurve sollte sich als Überlagerung dieser beiden Funktionen darstellen, was später in Abbildung \ref{MkVgl} auch gut zu erkennen ist. An dieser Stelle sei darauf hingewiesen, dass die Asynchronmaschine ein ähnliches Verhalten (\textit{KLOSS}´sche Formel \autocite{binder_bibel}) aufweist. 

Weiterhin soll der Kurzschlussversuch zur Einzelverlustbestimmung genutzt werden. Die gesamte Verlustleistung $P_\+{d}$ muss von der antreibenden Maschine aufgebracht werden, sodass sie direkt mit dem, an der Messwelle aufgenommen, Drehmoment $M_\+{k}$ bestimmt werden kann.
\begin{equation}
P_\+{d}= 2\pi M_\+{k}n
\label{eq:ku6}
\end{equation}
Es wird angenommen, dass sich die Verluste im Kurzschlussversuch aus den Stromwärmeverlusten $P_\+{Cu}$, den stromabhängigen Zusatzverlusten $P_\+{Cu,ad}$, den Reibungs- und Zusatzverlusten $P_\+{R+Z}$ und den Ummagnetisierungsverlusten $P_\+{Fe}$ (aufgrund der hohen Streuinduktivität $L_\+{s,$\sigma$}$ können diese nicht vernachlässigt werden) zusammensetzen:
\begin{equation}
P_\+{d}=P_\+{Cu}+P_\+{Cu,ad}+(P_\+{Fe})+P_\+{R+Z}
\label{eq:ku7}
\end{equation}
Die Stromwärmeverluste können bei bekanntem Strangwiderstand und Kurzschlussstrom berechnet werden:
\begin{equation}
P_\+{Cu}=3 R_\+{s} I_\+{s,k}^\+{2}
\label{eq:ku8}
\end{equation}
Die Reibungs- und Zusatzverluste können, wie später gezeigt wird, für jede Drehzahl aus dem Blindstromversuch entnommen werden. Um die Ummagnetisierungsverluste zu bestimmen, ist es notwendig den Betrag der Hauptfeldspannung $U_\+{h}$ zu kennen. Aus dem Ersatzschaltbild (Abbildung \ref{Ersatzschaltbild_kurz}) folgt:
\begin{equation}
U_\+{h}=\sqrt{R_\+{s}^\+{2}+(\omega_\+{s}L_\+{s,$\sigma$})^\+{2}}I_\+{s,k}
\label{eq:ku9}
\end{equation}
Für die Bestimmung der Hauptfeldspannung wird die Streuinduktivität $L_\+{s,$\sigma$}$ benötigt, welche aus der Geometrie der Maschine oder dem Bohrungsfeldversuch ermittelt werden kann. Sind aus dem Leerlaufversuch die Ummagnetisierungsverluste im unbelasteten Fall $P_\+{Fe,0}$ für jede Drehzahl bekannt, so können die Ummagnetisierungsverluste nach folgender Formel berechnet werden:
\begin{equation}
P_\+{Fe}(n)=P_\+{Fe,0}(n) \frac{U_\+{h}^\+{2}(n)}{U_\+{p}^\+{2}(n)}
\label{eq:ku10}
\end{equation}
Nun sind, bis auf die stromabhängigen Zusatzverluste, alle Einzelverluste bestimmt. Aus der Differenz der Gesamtverluste und der Einzelverluste können für jede Drehzahl die stromabhängigen Zusatzverluste errechnet werden:
\begin{equation}
P_\+{Cu,ad}(n)=P_\+{d}(n)-P_\+{Cu}(n)-P_\+{R+Z}(n)-P_\+{Fe}(n)
\label{eq:ku11}
\end{equation}
Werden die stromabhängigen Zusatzverluste durch die Stromwärmeverluste geteilt, ergibt sich der Faktor $k_\+{ad}$:
\begin{equation}
k_\+{ad}(n)= \frac{P_\+{Cu,ad}(n)}{P_\+{Cu}(n)}
\label{eq:ku12}
\end{equation}
Mit diesem Faktor können nun bei Kenntnis der Stromwärmeverluste für jede Drehzahl die stromabhängigen Zusatzverluste berechnet werden.
\section{Programm}
\begin{figure}[htb]
	\centering
  \includegraphics[width=1\textwidth]{auswertung_kurzschluss/kurzschluss_gui.png}
	\caption{Benutzeroberfläche zur Auswertung des Kurzschlussversuchs am Beispiel der Maschine EW2-A}
	\label{gui_kurz}
\end{figure}
Abbildung \ref{gui_kurz} zeigt die Benutzeroberfläche zur Auswertung des Kurzschlussversuchs. Die Benutzeroberfläche enthält zwei Achsensysteme, in denen die gemessenen Verläufe des Kurzschlussstroms (oberes Achsensystem) und des Kurzschlussmoments (unteres Axensystem) aufgetragen sind. Zur besseren Veranschaulichung sind beide Graphen bei gleicher Skalierung der Abszisse übereinander dargestellt. Außerdem können die beiden Achsensysteme über den Push Button \textit{In Figure öffnen} in zwei zusätzlichen Fenstern geöffnet werden. 

Weiterhin beinhaltet die Benutzeroberfläche ebenfalls einen Push Button \textit{Einstellungen}. Dort können die jeweiligen Spalten der Messdatei gewählt werden, in denen die, für diesen Versuch relevanten, Messwerte gespeichert sind.

Zur Berücksichtigung der temperaturabhängigen Größen steht eine Listbox im Panel \textit{Wahl des Temperatursensors} zu Verfügung. Es werden alle Temperatursensoren angezeigt, die bei Aufnahme der Messwerte vom Benutzer ausgewählt waren und es sollte derjenige gewählt werden, der die größte Übereinstimmung mit der Temperatur der Statorwicklung darstellt. Wurden mit einem Sensor aus technischen Gründen keine Messwerte aufgezeichnet, so erscheint eine Fehlermeldung. Zudem ist es möglich, mehrere Temperatursensoren auszuwählen. Dies führt dazu, dass zu jedem Messzeitpunkt die arithmetisch gemittelte Temperatur berücksichtigt wird.

Der Push-Button \textit{Auswertung} führt dazu, dass die Tabelle mit den gemessenen und den berechneten Größen geladen wird. Die erste Spalte enthält die Drehzahlen, die in der Messung angefahren wurden. Es sei darauf hingewiesen, dass jede Drehzahl aus Gründen der Übersichtlichkeit nur einmal in der Tabelle aufgeführt sind. Dazu wurden alle bei der Drehzahl gemessenen Größen arithmetisch gemittelt. Die zweite Spalte enthält die dazugehörige Frequenz, die intern aus Polpaarzahl und Drehzahl bestimmt wird. In der dritten und vierten Spalte befinden sich der gemessene Kurzschlussstrom und das gemessene Kurzschlussmoment. Die fünfte Spalte enthält den Betrag der Polradspannung $U_\+{p}$ , der aus der Drehzahl und dem aus dem Leerlaufversuch bekannten Faktor $k_\+{Up}$ berechnet wird. Die Polradspannung ist nötig, um die Induktivität $L_\+{d}$ nach Gleichung \ref{eq:ku3} zu bestimmen, welche in der sechsten Spalte eingetragen ist. Für diese Berechnung wird die Temperatur (gemäß der Auswahl des Temperatursensors) berücksichtigt und der Strangwiderstand nach Formel \ref{eq:widerstand} für jeden Messpunkt angepasst. In der siebsten bis neunten Spalte befinden sich berechnete Werte, die für die Gesamtverlustbestimmung in Abschnitt \ref{auswertung_verluste} benötigt werden: Der Betrag der Hauptfeldspannung $U_\+{h}$ nach Gleichung \ref{eq:ku9}, die gesamte Verlustleistung $P_\+{d}$ nach Gleichung \ref{eq:ku6} und die bei angepasstem Strangwiderstand bestimmten Stromwärmeverluste $P_\+{Cu}$ nach Gleichung \ref{eq:ku8}. In der letzten Spalte befinden sich Check-Boxen, die anzeigen, ob eine Zeile zur Mittelung berücksichtigt werden soll oder nicht. Diese Spalte ist folglich editierbar. Darüber hinaus können die Werte der Tabelle über den Push-Button \textit{Tabelle exportieren} in Excel exportiert werden.

Die letzte Spalte der Tabelle kann auch hier nicht nur manuell editiert werden. Zur automatischen Auswahl können die Funktionen des Panels \textit{Datenauswahl zur Mittelung} analog zur Auswertung des Leerlaufversuchs genutzt werden.

Das Panel \textit{Mitteln} dient zur arithmetischen Mittelung des Kurzschlussstroms $I_\+{s,k}$ und der Induktivität $L_\+{d}$. Dazu werden jeweils nur die, in der Spalte "`Berücksichtigen"' gewählten, Zeilen gemittelt. Weiterhin können $I_\+{s,k}$ und $L_\+{d}$ getrennt voneinander gemittelt werden, falls jeweils andere Zeilen berücksichtigt werden sollen.

Das Panel \textit{Fitting auf Drehmoment} besteht aus zwei Funktionalitäten. Die erste Funktionalität ist ein Algorithmus, der den Strangwiderstand $R_\+{s}$ und die Induktivität $L_\+{d}$ schätzt. Hierfür wird eine Funktion von Herrn Dipl.-Ing. Reis genutzt, die auf dieses Programm angepasst wurde. Die Funktion "`funcKurzschlussMomentFitting.m"' nutzt hauptsächlich die $\+{MATLAB}^\+{\textregistered}$-Funktion "`fminsearch(fun,x0)"', welche eine Optimierungsfunktion zum Finden des lokalen Minimums einer quadratischen Fehlerfunktion "`fun"' bei gegebenen Startwerten "`x0"' ist. Das Verfahren ist ein Downhill-Simplex-Verfahren nach \textit{Nelder-Mead}. In der Fehlerfunktion wird die Abweichung zwischen dem gemessenen Kurzschlussmoment $M_\+{k}$ und dem, nach Gleichung \ref{eq:ku5} berechneten, elektrischen Kurzschlussmoment $M_\+{e,k}$ gebildet und quadriert. Deshalb müssen an diese Funktion die gemessenen Drehmoment-Drehzahl-Paare, die Polpaarzahl und der Faktor $k_\+{Up}$ übergeben werden. Wird der Push Button \textit{Ausführen} betätigt, öffnet sich ein Dialogfeld, in dem die Startwerte eingegeben werden können. Nach Eingabe der Startwerte wird der Algorithmus gestartet und die Ergebnisse werden in den Textfeldern angezeigt. Die Ergebnisse oder manuell in die Textfelder eingetragene Werte können für die zweite Funktionalität genutzt werden. Nach Klicken auf den Push Button \textit{Anzeigen} wird in einem neu geöffneten Fenster der Vergleich zwischen den gemessenen Drehmomentwerten $M_\+{k}$ und den, nach Gleichung \ref{eq:ku5} mit den Parametern aus den Textfeldern berechneten, elektrischen Drehmomentwerten $M_\+{e,k}$ angezeigt. Die Anzeige ist beispielhaft auf Abbildung \ref{MkVgl} zu sehen. Gerade bei hohen Drehzahlen ist eine große Differenz zwischen beiden Drehmomentverläufen zu erkennen, was daran liegt, dass das gemessene Drehmoment um den Anteil, der durch Reibungs- und Zusatzverluste im Rotor $P_\+{R+Z}$ verursacht wird, größer ist. 
\begin{figure}[!h]
\begin{center}
\begin{tikzpicture}
\begin{axis}[
	x tick label style={/pgf/number format/.cd,%
         scaled y ticks = false,
					xtick={0, 500,  1000, 1500, 2000},
         set thousands separator={},
         fixed},
	ylabel=Drehmoment in \si{\newton\meter},	
	xlabel=Drehzahl in \si{\min^{-1}},
	%xunit=\metre\per\second, %This doesn't do what I want
	%y unit=\si{\metre\per\second}, %This does, but I don't want to type "\si{...}"
	%axis x line = bottom,
	%axis y line = left,
	xmin = 0, xmax = 2000,
	ymin = 0, ymax = 200,
	axis equal=false,
	grid=major,	
	legend pos=north west,
	height = 10cm,
	width = 1\textwidth]
	\legend{$M_\+{k,e}$ angepasst, $M_\+{k}$ gemessen}]
	%\addplot table {m_calc.txt};
	\addplot [style = solid,blue] table {auswertung_kurzschluss/m_calc.txt};
	\addplot [only marks, mark=star] table {auswertung_kurzschluss/m_mess.txt};
	%\addplot[style = solid] coordinates{(0,4) (1316,10.01)};
	% M = 4 + n*0.004568
\end{axis}
\end{tikzpicture}
\caption{Vergleich des berechneten und gemessenen Kurzschlussmoments am Beispiel der
Maschine EW2-A}
\label{MkVgl}
\end{center}
\end{figure}
%\begin{figure}[htb]
	%\centering
  %\includegraphics[width=1\textwidth]{auswertung_kurzschluss/M_k_vgl.png}
	%\caption{Vergleich des berechneten und gemessenen Kurzschlussmoments am Beispiel der Maschine EW-2A}
	%\label{MkVgl}
%\end{figure}
Abschließend können die Ergebnisse über den Push-Button \textit{In Projektdatei speichern} gespeichert werden. Standardmäßig werden der gemittelte Kurzschlussstrom $I_\+{k}$ und die gemittelte Induktivität $L_\+{d}$ in den vorgesehenen Parametern der Maschinen gesichert. Soll stattdessen der, aus der Optimierungsfunktion ermittelte, Wert für die Induktivität $L_\+{d}$ gespeichert werden, so ist die Checkbox \textit{gefittete Daten speichern} auszuwählen. Dies führt zusätzlich dazu, dass der Wert des Strangwiderstandes $R_\+{s}$ überschrieben wird. Außerdem werden für die Berechnung der Einzelverluste die Hauptfeldspannung $U_\+{h}$, die gesamten Verluste $P_\+{d}$ und die Stromwärmeverluste $P_\+{Cu}$ bei jeder gemessenen Drehzahl $n$ in der Variablen \textit{kurzschluss\_verluste} gespeichert.	
\section{Messergebnisse}

Tabelle \ref{tab:vgl_kurz} stellt die im Projektseminar gemessenen Werte denen aus \cite{deak} gegenüber. Es wird zum einen der Kurzschlussstrom $I_\+{k}$ und zum anderen die Längsinduktivität $L_\+{d}$ verglichen. Es ist zu erkennen, dass die Ergebnisse für $L_\+{d}$ sehr gut übereinstimmen. Die Abweichungen lassen sich durch Messabweichungen und die Mittelung begründen. Der Kurzschlussstrom $I_\+{k}$ stimmt bei der Maschine EW2-A wiederum sehr gut überein, wohingegen bei der Maschine I-C eine deutliche Abweichung festzustellen ist, was auf den reduzierten magnetischen Fluss als Folge des Demagnetisierungsversuchs aus \cite{deak} zurückzuführen ist, da unter Vernachlässigung des Strangwiderstandes und unter Berücksichtigung von Gleichung \ref{eq:le3} für den Kurzschlussstrom folgende Näherung gilt:
\begin{equation}
I_\+{k}\approx \frac{\Psi_\+{p}}{\sqrt{2} L_\+{d}}
\label{eq:ku13}
\end{equation}
\begin{table}[htbp]
  \centering
  \caption{Vergleich des Kurzschlussstroms $I_\+k$ und der Längsinduktivität $L_\+d$ der Maschinen EW2-A und I-C mit \cite{deak}}
    \begin{tabular}{l|cc|cc}
    \toprule
          & \multicolumn{2}{c|}{$I_\+{k}$ in \si{\ampere}} & \multicolumn{2}{c}{$L_\+{d}$ in \si{\milli\henry}} \\
    \midrule
    Maschine & EW2-A & I-C & EW2-A & I-C \\
    lt. \cite{deak}  & 76,8 & 71,3 & 2,78 & 2,45 \\
    Messung & 74,95 & 57,41 & 2,88 & 2,54 \\
    \bottomrule
    \end{tabular}%
  \label{tab:vgl_kurz}%
\end{table}%
Aus Gleichung \ref{eq:ku13} folgt unmittelbar, dass bei gleich bleibender Längsinduktivität und geringerem verkettetem magnetischem Fluss der Kurzschlussstrom geringer ist, was in Tabelle \ref{tab:vgl_kurz} deutlich wird. Darüber hinaus bietet es sich an, die Ergebnisse der Fitting-Funktion zu vergleichen. Tabelle \ref{tab:vgl_kurz2} enthält zum Vergleich den Widerstandswert aus der direkten Widerstandsmessung und die Längsinduktivität, die gemäß Gleichung \ref{eq:ku3} aus dem Kurzschlussversuch bestimmt wurde. Grundsätzlich ist zu sagen, dass sich die ermittelten Werte in der Nähe der Referenzwerte befinden. Die Abweichungen können sich mitunter dadurch ergeben, dass sich der Widerstandswert während der Messung temperaturabhängig verändert hat. Es ist als Erweiterung denkbar, eine Anpassung des Widerstandswerts in der Funktion vorzunehmen. Überdies hängen die Ergebnisse stark von den Messergebnissen ab, sodass Messfehler eine Ursache für abweichende Ergebnisse sein können. Wie in Abschnitt \ref{kurz_theo} bereits erläutert, ergeben sich gerade bei hohen Drehzahlen eine Abweichung zwischen berechnetem und gemessenen Kurzschlussmoment, sodass eventuell das Intervall, welches zum Vergleich herangezogen wird, angepasst werden muss. Abschließend sei noch einmal darauf hingewiesen, dass der zugrunde liegende Algorithmus lediglich ein lokales Minimum findet und die Lösung daher nicht eindeutig ist.
\begin{table}[htbp]
  \centering
  \caption{Vergleich des Strangwiderstands $R_\+s$ aus der Widerstandsmessung und der Längsinduktivität $L_\+d$ aus dem Kurzschlussversuch mit den Ergebnissen der Fitting-Funktion}
    \begin{tabular}{l|cc|cc}
    \toprule
          & \multicolumn{2}{c|}{$R_\+{s}$ in \si{\milli\ohm}} & \multicolumn{2}{c}{$L_\+{d}$ in \si{\milli\henry}} \\
    \midrule
    Maschine & EW2-A & I-C & EW2-A & I-C \\
    Fitting-Funktion  & 59.69 & 42.57 & 2.96 & 2.61 \\
    Widerstand/ Kurzschlussversuch & 55,46 & 43,12 & 2,88 & 2,54 \\
    \bottomrule
    \end{tabular}%
  \label{tab:vgl_kurz2}
\end{table}%

Abbildung \ref{kurz_Ik} stellt den Verlauf des Kurzschlussstroms über der Drehzahl für beide Maschinen dar und Abbildung \ref{kurz_Mk} bildet den Verlauf des Kurzschlussmoments über der Drehzahl ab. Die ausführlichen Auswertungstabellen sind im Anhang \ref{anhang_tabellen} aufgeführt. 
\begin{figure}[htb]
\begin{center}
\begin{tikzpicture}
\begin{axis}[
	x tick label style={/pgf/number format/.cd,%
          ytick={0, 20, 40, 60, 80},
					xtick={-3000, -2000, -1000, 0, 1000, 2000, 3000},
         set thousands separator={},
         fixed},
	ylabel=$I_\+{k}$ in \si{\ampere},	
	xlabel=$n$ in \si{\min^{-1}},
	%xunit=\metre\per\second, %This doesn't do what I want
	%y unit=\si{\metre\per\second}, %This does, but I don't want to type "\si{...}"
	%axis x line = bottom,
	%axis y line = left,
	xmin = -3000, xmax = 3000,
	ymin = 0, ymax = 90,
	axis equal=false,
	grid=major,	
	legend pos=south west,
	%height = 10cm,
	width = 0.8\textwidth,
	height = 8 cm]
	\legend{Maschine EW2-A, Maschine I-C}
	%\addplot table {m_calc.txt};
	\addplot [only marks, mark=triangle*] table {auswertung_kurzschluss/Ik2A.txt};
	\addplot [only marks, mark=*] table {auswertung_kurzschluss/IkIC.txt};
	%\addplot[style = solid] coordinates{(0,4) (1316,10.01)};
	% M = 4 + n*0.004568
\end{axis}
\end{tikzpicture}
\caption{Verlauf des Kurzschlussstroms der vermessenen Maschinen}
\label{kurz_Ik}
\end{center}
\end{figure}
\begin{figure}[htb]
\begin{center}
\begin{tikzpicture}
\begin{axis}[
	x tick label style={/pgf/number format/.cd,%
          %ytick={0, 20, 40, 60, 80},
					xtick={-3000, -2000, -1000, 0, 1000, 2000, 3000},
         set thousands separator={},
         fixed},
	ylabel=$M_\+{k}$ in \si{\newton\meter},	
	xlabel=$n$ in \si{\min^{-1}},
	%xunit=\metre\per\second, %This doesn't do what I want
	%y unit=\si{\metre\per\second}, %This does, but I don't want to type "\si{...}"
	%axis x line = bottom,
	%axis y line = left,
	xmin = -3000, xmax = 3000,
	ymin = 0, ymax = 200,
	axis equal=false,
	grid=major,	
	legend pos=north west,
	%height = 10cm,
	width = 0.8\textwidth,
	height = 8 cm]
	\legend{Maschine EW2-A, Maschine I-C}
	%\addplot table {m_calc.txt};
	\addplot [only marks, mark=triangle*] table {auswertung_kurzschluss/Mk2A.txt};
	\addplot [only marks, mark=*] table {auswertung_kurzschluss/MkIC.txt};
	%\addplot[style = solid] coordinates{(0,4) (1316,10.01)};
	% M = 4 + n*0.004568
\end{axis}
\end{tikzpicture}
\caption{Verlauf des Kurzschlussmoments der vermessenen Maschinen}
\label{kurz_Mk}
\end{center}
\end{figure}